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We investigate the optimal estimation of quantum expectation value of a physical observable, 
which minimizes a mean error with respect to general measure of deviation, when a finite number 
of copies of a pure state are prepared. If pure sates are uniformly distributed, the minimum value 
of mean error for any measure of deviation is achieved by projective measurement on each copy. 
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I. INTRODUCTION 

In the quantum information theory, it is one of im- 
portant issue to pick out the best way how to extract 
information from given quantum state. Many problems 
related to this have been investigated for a long time since 
Helstrom's and Holevo's works pl@]- In these problems, 
many copies of quantum state are needed, because we 
can get few information from one state, owing to sta- 
tistical properties of quantum mechanics. The measure- 
ments on these copies are classified into two types. One 
is separable measurement whose positive operator val- 
ued measure(POVM) is described by tensor products of 
POVM on each copy, the other is joint measurements 
whose POVM includes elements with entangled eigen- 
vectors. 

In many cases like quantum state estimation El 0, 0, 
sending some information by qubits 0, 0, and 
state identification 0,^2' ^ seems that the joint mea- 
surement is optimal. For example, in the paper 
Hayashi et al. discussed a problem of identifying a given 
pure state with one of two reference pure states, when no 
classical knowledge on the reference states is given, but 
a certain number of copies of them are presented. It was 
shown that the averaged success probability takes the 
maximum when we adopt joint measurement on tensor 
product of reference states and unknown state. However, 
Masahito Hayashi 0] investigated an estimation with 
minimum mean error with respect to general measure 
of deviation which satisfies appropriate conditions and 
showed that the minimum mean error over all measure- 
ment coincides with minimum mean error over separable 
measurement when the number of copies goes to infinity. 
For different situations [ill IIH fill . the similar result was 
obtained. 

In the paper D'Ariano et al. discussed the opti- 
mal estimation of an observable's expectation value for 
the finite copies of pure-state and showed that the mea- 
surement which minimizes the variance is separable one 
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if POVM is unbiased namely averaged estimator over the 
repeated measurement is equal to the expectation value. 
In the previous paper |l8[ , when N copies of pure state 
distributed uniformly in d-dimensional Hilbert space is 
presented, it was seen that the measurement which mini- 
mizes the averaged variance over all pure state is separa- 
ble measurement on each copy of N pure states. Thus for 
the estimation of observable's expectation value, it seems 
that separable measurement is optimal. In this paper, 
in order to make sure of this statement, we investigate 
whether minimum value of mean error with respect to 
any measure of deviation is achieved by separable mea- 
surement of each state or not. 



II. ESTIMATION WITH MINIMUM MEAN 
ERROR FOR GENERAL MEASURE OF 
DEVIATION 

In the previous paper [l8j , we determined the optimal 
way to estimate expectation value of a physical observ- 
able O, namely, POVM and estimator for each element 
of it which minimize the variance, when N copies of un- 
known pure state \4>) on d-dimensional Hilbert space are 
prepared. In the same situation we find the POVM {E a } 
and estimator uj a which minimize the mean error (A) 
defined by, 

< A > = (jy^\p m E a ]W(cJ a - Tr[n P ])^ , (I) 
where p is density matrix for pure sate \4>), 

bracket (• • ■ ) means the average for \<f>) which is dis- 
tributed uniformly over the d-dimensional Hilbert space 
and the function W(x) is a general measure of deviation. 
For simplicity, we assume that this function satisfies the 
conditions; 

(a) W{x) > (x + 0) and W(0) = 0, 

(b ) « > (x > 0), < o (x < 0) 



(c) 



d 2 W(x) 
dx' 2 



> 



2 



We will discuss the case where W(x) satisfies weaker con- 
ditions in the next section. 
Expanding the pure state \<j>), 



in the eigenvectors |n), (n = 1, 2, • • • d) with eigenvalue 
A„ of the observable fi, 

n\n) = A„|n) (n=l,2,-.- ,d), 

the density matrix p®-" becomes 

d d d d 

= ' ' ' ' ' ' c ni c mi ' ' ' c n N c m N 

ni = l n^ — lmi— 1 m^v = l 

x|m)(mi| ® ••• <g) |njv)(mjv|, 

- E E 

> o »i »J, > 

.1 + ••• + . „ = AT si+ ... + si = N 

x^-c5*(^)'i-.( C 3)'- 



E E 

(«!■• ' " ,*d) («!' ••• > 



x|ni)(mi| ® • ■ • ® |n-iv)(m.iv|- 



(2) 



Here, s n (n = 1,2, ••• ,d) is occupation number of 
the eigenstate \n) (n = 1,2, ■ ■ ■ , d) so the notation 
means that the summation is taken over 



»!.••■. «d > 
+ • ■ ■ + s d = N 



all non-negative integers s; (i = 1, 2, • • • , d) which sat- 
isfy the condition si + S2 + ■ ■ ■ + s c i = N and the notation 
means that the summation is taken over all 



states with same occupation number si (I = 1, 2, • • • , d). 

Substituting this expression for the density matrix p® N 
into the definition^) of (A) we want to minimize, we 
have 

< a > = |e/ n^nMEi^i 2 - 1 ) 

a •* n=l \n=l / 

x E E 



n ,...,s d > 
+ • ■ ■ + s d = N 



■ + s ' = N 



xc i ••• c /( c i) 1 •■•(cj d 



x E E 

{»!,-•• > "JV > {'"1 ■ ■ ■ ■ ■ m N } 

■■••d) (•!••••• «5) 

x|ni)(mi| ® • • • <g> |ra/v)(mjv| 



(3) 



where dc n dc* n stands for the integration over real part 
c„r and imaginary part c n \ of complex number c n and Z 



is normalisation factor, 



/e£ / d 

n dc„<«5 2 
n=\ \n=\ 

/d / d 

Y\ dc nR dc n iS ^2 



c„r - i 



,2 \ _ 1 27r 



Changing variables of integration c„ to £ n and y> n , 



c„ =C n e lip " (n = 1,2,- ■• ,d), 



the mean error (A) becomes 



^) - ^E E 



E 







a si, . . . ,s d > s > . . . , ,' > o 
>! + ■■■ +*d = » s' 1 + ...+5> d =N 

2tt / d \ 



\n=l 



E E 

{"1, ■ ■ ■ . " N > ■ ■ ■ . 

c«i. •••.«*) (.!,•••.. J) 

si+si + l f s d +s d + l 



x TrflmXrml 



C+Sd+ t^ 



|nAf)(mAf|S a ] 



n=l 



The contribution of the term where s n is different from 
s' n vanishes under the dip n (n — 1, 2, • • • , d) integration, 
and we have 



(A)= E 



su ... tad (u a )Tr[M Sl ... Sn E a }, (4) 



si, ... s d 



S1 + ■ ■ ■ + s d = N 



where the function w Slr .. jSd (w ) and operator A4 Sl ... s 
are defied by 



W Su ... , Sd (u a ) 



(27T) 



M sr .. s „ 



- 70 n=l \n=l / 

ne" +1 ^f^-E A ^«) ( 5 ) 

n=l \ n=l / 

E E . 



{n ll ■ ■ ■ ,n N } {m 1 , ■ ■ ■ , m N } 

Cn. - ■ ■ Cn.-" .= d ) 



x |m)(mi| ® ■ • ■ ® |njv)(mjv|- (6) 
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Clearly, the operator M. Sl ... Sn is non- negative and its 
rank is equal to one; 



Under this choice of POVM and estimators, the mean 
error (A) in equation becomes 



Ms 



\ 



E K> 



\ {»!,••• ,n N } 
\ ( = 1. ■ •• ,' d ) 



E 



(mi 



\n N ) 



<g> (mjv| 



From the conditions (a) ~ (c) for TU(:r), we can see 
that the first derivative of function w Slt ... , Sd (u> a ) becomes 
positive(negative) as x goes to +oo(— oo) and that the 
second derivative of it is positive. Thus there is only 
one value f2i™! n ? Sd where the first derivative of the func- 
tion w Sl .... . Sd (uJa) vanishes and the function iu Sl! ... jSd (w a ) 
takes the minimum. As the operator A4 Sl ... Sn is non- 
negative, we can get lower limit of mean error (A), 



(A) 



> 



E 

■ * » s d 
•• + «, 

E 



II* s 



+ ■ ■ ■ + s d = JV 



II* S 



n. ■■-.«*> o 

+ ■ ■ ■ + s d = JV 



x Tr[M sl ... Sn E a ], 
E w * 



Sd (u a )Tr[M Sl ... Sd E a ] 

^(min) « 



n. ■ ■ ■ ,*d > ° 

+ ■ ■ ■ + s d = JV 



x Tr M Sl ... Sn (j2 E ^ ' 

E/ (min) \ 

•1 'j>« 

n + ■ ■ ■ + s d = JV 

x Tr[M Sl ... s J, 



(7) 



(A) = E 



E 



•i,... ,»i > »l •i > 

.i + --- + *2, = w • ! +•■• + .,, = * 

x Wsi ,..., Sd (fiJ in) s ,)Tr[^ Sl ... s „7' si ,.., s ,] 

- E 

• ,. d > 

•1 + ' ' ■ +»d = N 

x w Su ..., Sd (n ( s ^\ Sd )Tv[M Sl ... s J, 
since V Sl ---, Sd and A4 Sl ... 5n satisfy the following relation 

Thus we can obtain the POVM {V Sl ..... Sd } and esti- 
mators Oi™! n ? )S(j which minimizes the mean error (A) for 
any measure of deviation . Eigenvectors of some elements 
of this POVM are entangled states of N copies (i.e any 
linear combination of states with same occupation num- 
ber) and the measurement for this POVM is joint one. 
However, in the light of the facts that elements of this 
POVM are linear combination of the projection opera- 
tors |ni)(ni| ® • • • l S> |?t-jv)(^jv| and that the mean error is 
linearly dependent on the element of POVM, we have an- 
other optimal measurement described by projection op- 
erators P ni ... nN (ni, • • • ,n N = 1,2, • • • ,d) 



\ni)(m\ ® ••• <8> \n N )(n N \, 



and estimator 171™™? Sd , where {si, • ■ • , Sd} is occupation 
number of the state after measurement, namely, the state 
where the i-th state is the eigenvector \rii), and this is 
separable measurement. 



Now, we prove that this lower limit is achieved if we 
choose V Sl - ,s d which is defined by 

"P sl -, Sd = E K)( n il ® • • • ® \ n N)(n N \, (8) 

{"!.■ ■ ,"JV> 

and which satisfies the condition for projective operators 
y , "Pat— ,s d = 1j 
V Sl -- ,s d T's l 1 --- ,s' d — 3s 1 s' 1 ' ' ' 5 aia ijV Sl — ,s d - 

as the POVM and fii™!"? Sd as the estimator for each 
element of this POVM. 



III. SUMMARY AND DISCUSSION 

In this paper, when iV copies of pure state which is 
distributed uniformly in rf-dimcnsional Hilbert space are 
prepared, we showed that the separable measurement is 
optimal in the sense that the mean error with respect to 
any measure of deviation which satisfies the condition (a) 
~ (c) takes the minimum. 

The POVM is independent of the choice of measure 
W{x) of deviation , although the estimator ils™! n .' ) iS<J is 
dependent on it. For the case with d — 2, N — 2 and 
f2 = o>, when we choose 

W(x) = er 2 sinh 2 -, (9) 
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as measure of deviation, the estimator to SltS2 is given by 

a ( {a 2 - 2a + 4) sinh \ - 2{a - 2) cosh § 
W2 '° = 4 ° g \ (a 2 +2(7 + 4) sinh | - 2{a + 2) cosh § 



w 0,2 





(j f ((7 2 + 2(7 + 4) sinh — 2(<r + 2) cosh — 
— los s S — 

4 ° 1 /_2 I- I /iNni-l, 2 Ol „or,l, 2 



(ct 2 -2(7 + 4) sinh f - 2(cr - 2) cosh % 



"^2,0 



which is dependent on the parameter a. 

The behavior of the estimator tt>2,o in the FIG.l shows 
that the limit of 0l>2,o as a approaches infinity is equal to 
i which is obtained from the result in the previous paper 
|18l | as is expected from the fact that the limit of measure 
of deviation W(x) as a tends to infinity is equal to x 2 . 




FIG. 1: The estimator LJ2.0 as a function of parameter a in 
measure© of deviation. When a goes to infinity, the estima- 
tor L02,o approaches ^ which is equal to the estimator for the 
case with d = 2, N = 2 and the observable o z in the previous 
paper [l^l . 

We assumed that the measure of deviation satisfies the 
conditions (a) ~ (c). These conditions are sufficient con- 
ditions for existence of unique value r^l™!"? Sd where the 
function w Slt ... . Sd (u> a ) defined by the equation (0 takes 
the minimum. Even if conditions (b) and (c) are replaced 
with weaker condition (b') 

(b') the function W(x) is continuous and monotoni- 
cally decreasing (increasing) function if x is neg- 
ative (positive), 

existence of minimum value of the function w Slr .. . Sd {ui a ) 
is shown as follows. Because the function w Sl! ... , Si (w a ) 
is continuous, in the finite region between the smallest 
eigenvalue A m ; n and the largest eigenvalue A max of ob- 
servable fi , there is, at least, one value f2'^ min - ) at which 
the function w Slt ... _ Sd (u> a ) takes the minimum. For the 
uj a which is smaller than the smallest eigenvalue A m i n of 
observable f2, because oj a — J2t=i i s negative; 



min < o, 



n=l 



the measure of deviation W(oj a — Yl n =i ^«£n) ^ s m ono- 
tonically decreasing function of variable tu a and the func- 
tion w sit ... , Sd (uj a ) is monotonically decreasing function 



of variable Lu a , too. Similarly, for the u a which is larger 
than the largest eigenvalue A max , w Slt ... , Sd (uj a ) is mono- 
tonically increasing function of variable uj a . Hence, at 
uj a = f2'^ mm '' this function becomes minimum over the 
interval (— oo,+oo). 

In the case with the condition (b') in stead of the con- 
ditions (b) and (c) , we may have some points at which 



the function 



,(wa) take the minimum and can- 



not uniquely determine estimator corresponding to each 
element of projection operator. 

We investigated the case where pure state is dis- 
tributed uniformly. Because of this assumption, inte- 
grand in the equation Q except for the factor Y\fj c ni c* n . 
is independent of complex argument Lp.i of coefficient c n 
and we obtained the equation (0J . More generally, for the 
case where pure states are distributed by the probability 
p(|ci|, ■ ■ • , \cn\), the mean error (A') becomes 

< a '> = e / ri^c^ncic n i 2 -i) 

a ^ n=l \n=l ) 

d d d d 

x P (\ Cl \,... ,\ CN \)J2--- E E ••• E 

ni — 1 njv — 1 mi — 1 mjv = l 

N 

x II Cn ' C ™. Tr [l ni ^ mi l ® l"iv>(miv|^ ] 

id 

d 

x IFK-^A„|c„| 2 ), 

n=l 

and we obtain equation Q with w' si ... s (oj a ) 

poo d / d \ 

Jo n=l Vn=l / 

d 

x p(ei,-,&on# n+1 

n=l 

instead of w Sl .... . Sd (uj a ). In the same manner as in the 
previous section, it is shown that the separable measure- 
ment on each copy is optimal. For example, we consider 
measurement of observable a z on 2-dimensional Hilbcrt 
space. Using polar coordinate (9, <p) on Bloch sphere, any 
pure state \4>) is described in the form, 

!</»)= cos ||0> sin 

up to phase factor, where |0) and |1) are eigenvectors of 
@ z 

(7 Z |Q) = +1|0), a z \\) = -l|l>. 

If the probability density for distribution of pure state is 
independent of azimuthal angle namely, if the distribu- 
tion of pure state is symmetric under the rotation around 
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z axis, separable measurement on each copy makes the ments. 
mean error minimum and is one of the optimal measure- 
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